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We solve the long-standing problem of making n perfect clones from m copies of one of two 
known pure states with minimum failure probability in the general case where the known states have 
arbitrary a priori probabilities. The solution emerges from a geometric formulation of the problem. 
This formulation also reveals a deeper connection between cloning and state discrimination. The 
convergence of cloning to state discrimination as the number of clones goes to infinity exhibits a 
phenomenon analogous to a second order symmetry breaking phase transition. 
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It is impossible to always clone quantum states per¬ 
fectly mm- This leads to advantages for quantum sys¬ 
tems over classical ones in quantum based communica¬ 
tions protocols of practical relevance. A common exam¬ 
ple is provable security in quantum key distribution, with 
recent works showing more applications [ 311 ]. Develop¬ 
ments in cloning, including deterministic but approxi¬ 
mate cloning I3{ini and probabilistically perfect cloning 
provide anchors for better understanding quan¬ 
tum theory as a whole, such as the relationship between 
the no-cloning and no-signaling theorems mi, and fun¬ 
damental limits on quantum measurements [T5HT3 . In 
particular, cloning’s relationship with the fundamental 
limits of state discrimination will be central to this Let¬ 
ter. For reviews citing recent developments, applications 
and experiments related to cloning see [M HO] and to 
state discrimination see m- 

When knowledge of the state preparation is available, 
perfect cloning is probabilistically possible. With the first 
result in probabilistic cloning, Duan and Guo m con¬ 
sidered the problem of producing perfect clones of lin¬ 
early independent pure states and focused on the two 
state case. They found the maximum average success 
rate when both states are equally likely, the case of equal 
priors, and set this success probability as lower bound 
for arbitrary prior probabilities. While other work has 
been done on this problem, there has until now been no 
general solution. In this Letter we obtain the general ana¬ 
lytic solution and, with its help, examine the relationship 
of cloning and state discrimination. 

There are a number of reasons why one might want to 
solve the general problem with arbitrary priors, (i) The 
solution to the equal prior problem is obtained using 
only symmetry arguments, with no need for optimiza¬ 
tion. (ii) A general solution would check the robust¬ 
ness of the case of equal priors against variations of the 
prior probabilities around 1/2. This gives control over 
errors that are unavoidable for any physical realization. 


(iii) One could consider a discrimination protocol con¬ 
sisting of optimal cloning followed by optimal Unambigu¬ 
ous Discrimination (UD) of the produced clones. This is 
a particular case of a “prepare and measure protocol” 
which we will call “discrimination by cloning.” Surpris¬ 
ingly, this is optimal for equal priors and for any number 
of clones, n, as shown later. This suggests that the case 
of equal priors is very special and can provide a decep¬ 
tive view of cloning, (iv) In the limit of infinitely many 
clones, the optimal strategy prepares the clones according 
to the outcomes of UD of the input states. This is a par¬ 
ticular case of a “measure and prepare” protocol, which 
we will call “cloning by discrimination.” Since the UD 
measurement varies over the range of prior probabilities 
(a 3-outcome generalized measurement vs. a 2-outcome 
projective measurement, |21jl. this hints at a similar sit¬ 
uation for optimal cloning that can only be revealed by 
solving the general problem. 

Our solution shows that discrimination by cloning as 
outlined in (iii) is sub-optimal for unequal prior proba¬ 
bilities (unless one state is never sent). This indicates 
that the equal prior case is not representative of state 
dependent cloning. Additionally, contrary to the sugges¬ 
tions in (iv) above, our solution leads to a failure prob¬ 
ability that is a smooth function of the priors. How¬ 
ever, the strategy converges to cloning by discrimination 
as n —>■ oo, implying a discontinuous second derivative 
and revealing a phenomenon similar to a second order 
symmetry breaking phase transition. 

We envision a state dependent probabilistic doner as a 
machine with an input port, an output port and two flags 
that herald the success or failure of cloning. The input 
IV’)”) = IV"!)®™! i = 1,2 (to identical copies of either 
\ipi) or |'!/’ 2 )) is fed through the input port for processing. 
In case of success, n perfect clones = It/i)®" are 
delivered through the output port with probability pi, 
conditioned on the input state being ItpT')- Otherwise, 
the output is in a generic failure state, with a failure 
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probability qi = 1 — pi. 

For cloning, optimality is usually addressed from a 
Bayesian viewpoint that assumes the states to be cloned 
are prepared with some prior probabilities rji and 772 , 
771+772 = 1- Then a natural cost function for the proba¬ 
bilistic cloning machine is the average failure probability, 

Q = Viqi + r] 2 q 2 - (1) 

The aim is to find the optimal doner that minimizes the 
cost function Q, and yields the minimum average failure 
probability Qmin for arbitrary priors 771 and 772 . 

In our formulation, similar to that in m, the Hilbert 
space IH®"* of the original m copies is supplemented by 
an ancillary space 0 Jfj? that accommodates 

the additional n — m clones and the success/failure flags. 
Next, we introduce a unitary transformation U via 

C/|7/’r)|0) = v^lO|a.)+v^J‘i>”)K), 7 = 1,2. (2) 

Here the ancillas are initialized in a reference state |0). 
The states of the flag associated with successful cloning, 
joi), are orthogonal to the state associated with failure, 
|ao), and jd)") is a generic failure state in IK®", the same 
for both inputs for optimality. 

We note that, although |ai) = |a 2 ) for optimal 
cloning, we consider a more general scenario allowing 
for |ai) ^ [ 02 ) to include “cloning by discrimination.” 
In this protocol we employ UD to identify the input 
state and then prepare clones of the identified state. For 
UD the success flag states must be fully distinguishable, 
so (Q!i|a 2 ) = 0. Further, an even more general sce¬ 
nario could allow for two different failure states, |<I>") 
(7 = 1,2), in Eqs. ([^. This, however, would necessarily 
be sub-optimal since we could probabilistically determine 
whether we received or \il)2) by applying UD to the 
states |$"). In case of success, we could prepare n copies 
of the state, increasing the overall success rate of the 
cloning strategy. 

To proceed, we now take the inner product of equa¬ 
tion i with itself in ([^, yielding that the probabilities 
are normalized, Pi + qt = 1. Taking the inner product 
of equation 7 = 1 with 7 = 2 in (H) yields the main con¬ 
straint, 

s'" = VP 1 P 2 s"a + VK'?2, (3) 

which is a consequence of the unitarity of U. Here we 
introduced the notation s = {'4’i\'4’2) and a = (ai|a2)- 
The overlaps s and a can be chosen to be real without 
any loss of generality, so 0 < s, a < 1. Clearly, a = 1 for 
optimal cloning, while a = 0 for cloning by discrimina¬ 
tion. If Eq. ^ is satisfied, it is not hard to prove that U 
has a unitary extension on the whole space. 

Before developing the analytical theory of optimizing 
(minimizing) Q, we present a complete geometric picture, 
similar in spirit to that in [^, that visually solves the 


optimization problem and serves as guide for the subse¬ 
quent calculations. For the geometrization the following 
features of ([^ turn out to be important, (a) For fixed s, 
n and m ([^ defines a class of smooth curves on the unit 
square 0 < < 7 ^ < 1 (e.g., solid, dashed or dotted curves 
in Fig. [^. (b) All these curves meet at their endpoints, 
(1, s^"") and (s^™, 1). (c) At the endpoints the curves be¬ 
come tangent to the vertical and horizontal lines <71 = 1 
and (72 = 1, respectively, provided a ^ 0. (d) For a = 0 
the curve defined by Eq. (|^ is the hyperbola ( 71(72 = 
(dashed line in Fig. [^. (e) The curve corresponding to 
a particular value of a and the segments joining its end 
points with the vertex (1,1) form the boundary of the 
set Sa (any of the gray regions in Fig.[^, where 

Sa = {((?i, (72) : ^/PiP2 + V'?i92 - s'" > 0 }. ( 4 ) 

The sets satisfy Sa C Sa' if a < a', (f) The sets Sa are 
convex if o; > 0. In particular Si is convex. 

These considerations lead to the emergence of a geo¬ 
metrical picture of the optimization problem which we 
display in Fig. 



FIG. 1: Unitarity curves, 52 vs. qi, from Eq. <§ and the 
associated sets Sa from Eq. Q for q = 1 (solid/light gray), 
a = 0.5 (dotted/medium gray), and o, = 0 (dashed/dark 
gray). The figure also shows the optimal straight segment 
(3 = 771 ( 71 - 1 - 772(72 and its normal vector ( 771 , 772 ). For all 
plots s = 0.5, 7Ti = 1, and n = 2. 

Eq. Q defines a straight segment on the unit square 
0 < (7i < 1 with a normal vector in the first quadrant 
parallel to (771,772). For fixed priors, the average failure 
probability Q is proportional to the distance from this 
segment to the origin (0,0). Since ^i is convex and the 
stretch of its boundary given by Eq. ^ with o; = 1 is 
smooth, a unique point ((71,52) of tangency with the seg¬ 
ment 0 exists for any value of the priors and finite n. 
It gives Qmin and defines the optimal cloning strategy. 

We note that the inclusion hierarchy of the sets Sa pro¬ 
vides a simple geometrical proof that a = 1, i.e., |ai) = 
\a 2 ), is the optimal choice for cloning. For “cloning by 
discrimination,” on the other hand, (ai|a 2 ) = a = 0. 
From points (d) and (e) above, it follows that for any 
finite n this protocol is strictly suboptimal, i.e., (3min < 
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Qud , noticing that the failure rate of cloning by discrim¬ 
ination is that of UD, Qud- This is intuitively obvious 
since in cloning one is asking for less than in UD; the 
identity of the input states is not revealed for any fi¬ 
nite n. However, optimal cloning and UD become one 
and the same in the limit n —>■ oo, when s" —>■ 0 and 
the curve ([^ collapses to the hyperbola qiq 2 = as 
it does for a = 0. We return to this point below. 

A more quantitative analysis requires finding a conve¬ 
nient parametrization of the constraint <§■ To this end, 
we set a = 1 and write = sinOi for 0 < 9i < 7r/2. By 
further introducing the variables x = cos( 0 i-|-d 2 ) and y = 
cos(0i — 62 ) we manage to linearize the curve (|^, which 
now is the straight segment 2s™ = (1 + s'^)y — (1 — s^)x 
with \x\ < y < 1. Its guiding vector is readily seen to 
be (1 -I- s", 1 — s"), so the segment’s parametric equation 
can be written as 

l-(l + s")t l-(l-s")t 

X =-, y =-, (5) 

gn—m ' gTi—m ’ ^ ' 

where we have rescaled the guiding vector so that Eqs. Q 
and below are simplest. Because of the symmetry of 
this procedure, the parameters x and y are invariant un¬ 
der <71 o (72 (equivalently, under 0i -o- 62 ). Thus, the two 
mirror halves of the curve (§ under this transformation 
are mapped onto the same straight line <§■ By express¬ 
ing qi as a function of t only half of the original curve is 
recovered. The other half is trivially obtained by apply¬ 
ing < 7 i o 52 - After putting the various pieces together 
one can easily get rid of the trigonometric functions and 
express Eq. (§ in parametric form as 

l-xy- (-l)Vl - xVl - 2/^ ■ 10 

qt = -^-, 1 = 1 , 2 . (6) 

Fig. i shows examples of the unitary curve @. 

(a) (b) 




FIG. 2: Unitarity curves, 52 vs. qi, from Eq. for different 
values of s and for (a.) m — 1, n = 2 and (b) m = 1, n = 5. 
The curves are symmetric under mirror reflection along the 
dotted line qi = ( 72 , i.e., under the transformation qi e-)- 52 - 
The dashed lines in (b) are the hyperbolae qiq 2 = 


for s close to 1 the hyperbolas remain closer to the 
vertex (1,1), but retain the same end points. In the 
limit n —00 all these curves become hyperbolic. 

We now return to finding the minimum, Qmin) of the 
average failure probability Q. Despite its apparent sim¬ 
plicity, this involves solving a high-order equation with¬ 
out a simple form. Instead, we will derive a parametric 
equation for Qniinivi)- Along with the complete descrip¬ 
tion of the unitary curve (§, this provides a complete 
solution of the problem in parametric form. 

With no loss of generality we may assume rji < 772 or, 
equivalently, 0 < 771 < 1/2. Then the slope of the straight 
line Q, — 771 / 772 , satisfies — 1 < — 771/772 < 0. Hence, it 
can only become tangent to the lower half of the unitarity 
curve (§ (see Fig. [^. Increasing qi, the slope of this 
lower half increases monotonically from —1 at qi = 52 , 
to 0 before we reach the line <71 = 1 (assuming n is finite). 
This follows from the properties (a)-(f) above and can be 
checked using Eq. Q. The values of t at which the slope 
is —1 and 0 are, respectively. 


t-i — 


1 - 


1 — s’^ 


1 _ g2(n-m) 

^0 = 1 - S2n • 


For any point {qi{t), q 2 {t)) with t G [t_i,to] there is a 
line Q = rjiqi + 772(72 that is tangent to it, starting with 
Vi = V 2 = 1/2 for t = t_i up to 771 = 0, 772 = 1 for t = to- 
This observation enables us to derive the desired para¬ 
metric expression for the optimality curve (5min(??i) as 
follows. For a given t in the range above, a neces¬ 
sary condition for tangency is r]iq[ + 772(72 = 0 , where 
( 7 ' = dqi/dt. In this equation we can solve for 771 (or 
772 ) using that rji + r ]2 = 1- By substituting qi and (72 
in Eq. Q with ([^ we enforce contact with the unitarity 
curve and obtain the expression of Qmin- The final result 
can be cast as: 


m 


92 

92 - 9i ’ 


Qn 


9^1 - 9i92 
92 -9i 


, t-i < t < to, (8) 


where t-i, to and qi are given in Eqs. Q and ([^. The 
expressions for the derivatives g' are 


9i = 


Vqtjl - di) 

gn—m 


l-bs" 
l/l — x"^ 


{-ly 


1 - s” 
\/i - 


(9) 


Fig. 1^ shows plots of the curves (5min(’7i) for tti = 1 
input copies and (a) n = 2 or (b) tt, = 5 clones, as in the 
previous figure. We see that Qmin is an increasing func¬ 
tion of 771 in the given range [0,1/2]. The values of Qmin 
at the end points of this range follow by substituting to 
and t_i, Eq. into Eq. They are given by 

2m _ 2n m _ n 

Qo = 92(to) = ^3^, = (10) 


For n > 2 the curves closely approximate qiq 2 = 
(dashed lines) for small and moderate values of s, while 


where Qmin = Q-i holds for equal priors and Qmin = Qo 
for 771 —7- 0 (i.e., 772 —)■ 1). The dashed lines in Fig. i(b) 
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FIG. 3 : Minimum cloning failure probability Qmin vs. rji 
(solid lines) and UD failure probability Quo vs. 771 (dashed 
lines) for the same values of m, n and s used in the previous 
figure. 


depict the well known optimal UD solution 


1 


o2m 


m 


S^'"772, 


0 < 77 i < 


1 + 


( 11 ) 


It is apparent from these plots that the optimal cloning 
protocol performs strictly better than cloning by discrim¬ 
ination, as was geometricaly proved in Figs. and 
However, the difference in performance decreases with in¬ 
creasing number of clones. In Fig. (b), for only n = 5 , 
a difference is hardly noticeable for s < 0 . 5 . For s > 0.5 
the convergence is slower but in the limit n —>■ 00 there 
is perfect agreement for any s < I. 

The convergence of the optimal cloning failure proba¬ 
bility, Qmin, to that of cloning by discrimination, Qud in 
Eq. ^ follows from our geometrical approach. Recall 
that in the limit n —> 00 (or equivalently a —>■ 0) the 
right hand side of Eq. ^ describes hyperbolas that we 
can write as q2 = (dashed lines in Figs. and [^. 

Their slopes are in the range [— 1 , — s^"]. A unique point 
of tangency with the line ([^ can only exists if the slope 
of this line, —771/772, is within this same range. This gives 
the 771 interval in the first line of (11), and one can easily 


obtain the corresponding expression for Quo- If the slope 
of the line ([^ is outside the range, tangency is not pos¬ 
sible, and the optimal line merely touches the end points 
of the hyperbolas, so the expression of Qud becomes the 


second line of Eq. ( 11 ). In geometrical terms, the straight 
line 0 pivots on the end points as we vary 771. Fur¬ 
thermore, note that for the second line in (11) we have 


Pi = 1 — <71 = 0, which leads to a 2-outcome projec¬ 
tive measurement, as only one success flag state (|a2)) is 
needed in Eqs. ([^. 

Interestingly, in this limit a phenomenon analogous to 
a second order phase transition takes place. Our geo¬ 
metrical approach shows that the average failure proba¬ 
bility Qmin( 77 i) is an infinitely differentiable function of 771 


for finite n. However, as n goes to infinity (or a —>■ 0 ) 
the limiting function (5ud(77i) has a discontinuous second 
derivative. Moreover, the symmetry gi o (72 breaks in 
the “phase” corresponding to the second line in Eq. ( 11 1. 
A similar phenomenon arises in UD of more than two 
pure states [ 55 ]. 

It has been argued above that cloning by discrimina¬ 
tion is strictly suboptimal (unless n —>■ 00). One could 
likewise wonder if discrimination by cloning can be op¬ 
timal. On heuristic grounds, one should not expect this 
to be so, as cloning involves a measurement and some 
information can be drawn from the observed outcome. 
However, the equal-prior and the 771 —)■ 0 cases provide 
remarkable exceptions. For both we may write the total 
failure rate as Qc + (1 — Qc)Qud, where C stands for 
cloning. For 771 = 772 = 1/2, Eq. ([^ implies Qc = Q-i, 
in which case the produced 7 i-clone states are equally 
likely. The UD of these states fails with probability s", 
as follows from Eq. (11) applied to n copies. The total 
failure rate is then s’", which is the optimal UD failure 
rate of the original input states, Eq. ( 11 ). If 771 —)■ 0 then 
only \1jj2) is produced with non-vanishing probability and 
Qc = Qo- Failure in the second step (UD) is given by 
the bottom line in Eq. 0 applied to n copies. The total 
failure rate is s^"", also achieving optimality. 

Using our main result in Eqs. (|^, and 0 one can 
check that these are the only cases where discrimination 
by cloning is optimal. These are also the only cases where 
no information gain can be drawn from the cloning mea¬ 
surement. This hints at how special these cases are and 
justifies the need of the derived solution for arbitrary pri¬ 
ors to have a full account of two-state cloning. 

In summary, we have provided the general solution to 
the long-standing probabilistically perfect cloning prob¬ 
lem for two states. The unequal prior case (771 7^ 772) 
uncovers remarkable phenomena that the very special 
equal priors case was unable to reveal. In particular, 
the convergence of cloning to unambiguous discrimina¬ 
tion as the number of clones becomes very large involves 
a discontinuity in the second derivative of the failure 
rate Qmin( 77 i), corresponding to a second order symme¬ 
try breaking phase transition. Our geometrical approach 
proved very helpful for both visualizing what the solu¬ 
tion looks like qualitatively and for deriving the analyt¬ 
ical solution. The same approach can also be applied 
to cloning of three or more states, deterministically ap¬ 
proximate cloning, and other optimization problems that 
involve highly non-linear constraints. 
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